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Abstract. We establish a fundamental connection between the tropical RSK correspondence and GL{N, R)- 
Whittaker functions, analogous to the well known relationship between the RSK correspondence and Schur 
functions. This gives rise to a natural family of measures associated with GL ( Af , R)- Whittaker functions 
which are the analogues in this setting of the Schur measures on integer partitions. The corresponding ana- 
logue of the Cauchy-Littlewood identity can be seen as a generalisation of an integral identity for GL{N,M.)- 
Whittaker functions due to Bump and Stade. As an application, we obtain an explicit integral formula 
for the Laplace transform of the law of the partition function associated with a one-dimensional directed 
polymer model with log-gamma weights recently introduced by one of the authors (TS). 



1. Introduction 

The Robinson-Schensted-Knuth (RSK) correspondence is a combinatorial mapping which plays a funda- 
mental role in the theory of Young tableaux, symmetric functions and representation theory. In particular, 
it provides a direct combinatorial proof of the Cauchy-Littlewood identity 



A ij 

where the sum is over integer partitions and s\ denotes the Schur function associated with the partition A. 
The Schur fmiction s\{x) is a symmetric function in the variables x = {xi,X2-, ■ • •) defined by 

T 

where the sum is over semistandard tableaux T with shape A and where /Xi is the number 

of i's in T. For more background on symmetric functions we refer the reader to [38]. 

The RSK mapping is defined by a combinatorial algorithm which associates to each matrix M = {rriij} 
with non-negative integer entries a pair (P, Q) of semi-standard tableaux with the same shape. By a version 
of Greene's theorem [30, 34], it can also be defined via expressions in the (max, -|-) semi-ring. Starting 
with the latter definition, A.N. Kirillov [34] considered a geometric lifting of the RSK mapping, obtained 
by replacing (max, by (-I-, x) in these expressions, which he called the tropical RSK correspondence. 
This has been studied further from a dynamical point of view by Noumi and Yamada [41] and related to 
Dodgson's condensation method for computing determinants in [19]. See also [7] for related work. 

It is natural to ask if there is an analogue of the Cauchy-Littlewood identity for the tropical RSK cor- 
respondence. A priori it is not at all clear that this should exist: something quite remarkable needs to 
happen in order to see the product structures on both sides which characterise this identity. In this paper 
we show that in fact these product structures do appear in the correct formulation and, moreover, the role 
of the Schur functions is now played by GL(iV, M)-Whittaker functions. In a particular case (corresponding 
to square matrices) the analogue of the Cauchy-Littlewood identity is in fact a well known integral identity 
for G'L(A^, M)-Whittaker functions due to Bump and Stade [15, 58, 27]. 

Before stating our main results, we will first explain the origin of these product structures in the context 
of the RSK correspondence and its role in the proof of the Cauchy-Littlewood identity. The RSK mapping 
associates to each matrix M = {mij} with non- negative integer entries a pair {P,Q) of semi-standard 
tableaux with the same shape. If M is an n x matrix then P has entries from {1,2,..., and Q has 
entries from {1, 2, . . . , n}. Moreover, the vector (i?i, . . . , i?„) of row sums of M is the type of P, that is, i?i 
is the number of I's in the semistandard tableau P, i?2 is the number of 2's, and so on. Similarly, the vector 
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of column sums (Ci, . . . , Cjv) is the type of Q. Note that this reflects the well known symmetry property of 
RSK, namely that if M ^ (P, Q) then M* (Q, P). It follows that 

and summing both sides gives the Cauchy-Littlewood identity: on the left we sum over all matrices with 
non-negative integer entries and on the right we sum over all pairs of semistandard tableaux P and Q with 
the same shape. 

Another way of interpreting the above argument is as follows. Let pi, . . . ,p„ and qi, . . . , be a collection 
of positive numbers such that < piQj < 1 for all i and j, and consider the probability measure on n x 
matrices defined by 

From the above discussion, the push-forward of this probability measure onto the shape of the tableaux 
obtained under the RSK mapping is given by 

(1-1) Pm)=l[(^-MMp)sx{q). 

Now, the Cauchy-Littlewood identity is essentially equivalent to the fact that P is a probability measure on 
the set of integer partitions. Such probability measures are known as Schur measures [4'.)]. 

We now turn to the tropical RSK correspondence. The input is an n x iV matrix X — {xy} with strictly 
positive real entries and, supposing here for convenience that n > N, the analogue of the P-tableau is 
a triangular array {zk,e)i<e<k<N of non-negative real numbers. The 'shape' of P is the vector zjv,-- Let 
9i, . . . ,9n and 9i, . . . ,6^ be a collection of real numbers satisfying 6i + 9j > for all i and j, and consider 
the product measure on input matrices X defined by 

where denotes the distribution of the inverse of a Gamma random variable 

The main result of this paper (Theorem 3.8) is that the push-forward of the probability measure fi under 
the tropical RSK mapping onto the 'shape' {zN,e)i<e<N is given by 

^ r(^m + Xi) -Q dyi 

where the functions 4'i/(y) are Gi(A^, ]R)-Whittaker functions (defined in Section 3.2 below) and 

^A^(A)=(^^nr(A.-A.)-. 

The probability measure /i^ is the analogue of the Schur measure in this setting and shall be referred to as 
a Whittaker measure; the fact that it integrates to one is the analogue of the Cauchy-Littlewood identity. 
In the special case iV = n, we obtain the following simplification: 

N N N , 

(1-3) f^%{dy) = n l[^{^^ + Lr'e~y-'^e{y)^g{y)Y[^. 

m=li=l i=l 

In this case, the analogue of the Cauchy-Littlewood identity reduces to 

N , N N 



(1.2) f^^{dy) = My)f dXs^{X)^^,{y)l[l[^^^^^ Yl 



/ e-^«'vl;,(y)vl,.(y)j]^= [][]r(0, + ^„), 

"'*+ i=l ^' m=li=l 



which can be seen to be equivalent to a Whittaker integral identity due to Bump and Stade [15, 58, 27]. We 
remark that, while this gives a new interpretation of the Bump-Stade Whittaker integral identity, we do not 
give a new proof of this identity in the present paper; in fact, we use it to obtain the simplification (1.3). 

We prove Theorem 3.8 by considering a dynamical version of the tropical RSK construction due to Noumi 
and Yamada [41], allowing n to increase as we successively add rows to the input matrix X. The image 
triangular array z{n) evolves as a Markov process in discrete time n subject to a particular entrance law for 
n < N. First we prove that the shape ZN^.{n) evolves marginally as a Markov process in its own filtration. 
The Markov property of ZM,-{n) relies (via the theory of Markov functions) on an algebraic intertwining 
relation between Markov kernels for z{n) and ZN^.{n) as well as on a limiting argument which shows that the 
result holds for the particular entrance law for z{n) dictated by the tropical RSK correspondence. Secondly 
we prove that this Markov process can be diagonalized in terms of G'L(A^, R)- Whittaker functions. This 
yields the formula (1.2) for the fixed n probability distribution of the shape. 

Our results have an important application to the study of directed polymers, analogous to the role of the 
RSK correspondence in the study of longest increasing subsequences and last passage percolation [2, 3, 31, 
32, 37, 60]. According to the definition of tropical RSK, we can write ZN,i{n) = X^tt j)g7r -^u where the 
summation is over all 'up/right' lattice paths in Z-^ from (1, 1) to {n,N). Under the measure fi, this random 
variable can be interpreted as a partition function for a directed polymer in a random environment given by 
the weights Xij. For a particular (homogeneous) choice of the parameters 9i + Oj =7 for all i and j, this 
model was introduced and studied in the paper ["i(i]. In particular, setting Zn = z„^i(n), it was shown in 
[56] that the free energy is given explicitly by 

(1.4) lim ilogZ„ = -2*(7/2) 

almost surely, where *I'(.t) = [logr]'(a;) is the digamma function, and moreover 

n K\ r varlogZ„ 

(1.5) limsup YR < 00. 

These asymptotic results were obtained via the observation that for this polymer model there is an analogue 
of the output theorem (or Burke property) for the single server queue. This observation, combined with 
recent progress on a related model in [42], provided the inspiration for the present work. In fact, the model 
considered in [42] is a degeneration of the one we presently consider (c.f. Section 4.2). 

Using our main result, we obtain the following explicit formula for the Laplace transform of the law of 
the partition function ZAr_i(n) under the measure /i (this statement is also contained in Theorem 3.8): 

(1.6) ^H = / dA.s~.(«-^-) n r(A.-^^,) nn^iS^^-A^W' 

where the poles of the functions r(Ai — dj) and r(Ai + 0,„) are not encountered as we may assume without 
loss of generality that Qm > for all m and 9j < for all j. 

Our formula (1.6) has recently been applied in [11] to prove the following asymptotic result: there exists 
7* > such that the inverse-gamma weight polymer free energy with parameter 7 £ (0,7*) has limiting 
fluctuation distribution given by 

where f-y = -2^(7/2), g-^ = -2*"(7/2) and i^cuE is the CUE Tracy-Widom distribution function. The 
restriction on the parameter 7 is present for purely technical reasons. 

There have been many recent developments on the Kardar-Parisi-Zhang (KPZ) equation [17, 24, 52] and 
various discretisations of this equation which have an underlying algebraic structure. The latter include the 
Whittaker measures (and processes) introduced in the present paper (and a particular family of degenerations 
of these which were introduced earlier in [42]). These in turn can be seen as degenerations of the Macdonald 
measures and processes introduced in [10] which gives access to the theory of Macdonald polynomials and 
incorporates a wide range of other interesting degenerations including g-deformations of the Whittaker 
measures and processes discussed here. A feature which is important about the Whittaker measures we 
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consider here is the fundamental hnk to the tropical RSK correspondence. In particular, our main result 
is not a degeneration of any of the main results in the more recent paper [10]. In another direction, it is 
possible to define an analogue of the tropical RSK mapping directly in the continuum setting of the KPZ 
equation, with input given by a two-dimensional strip of space-time white noise [46]. Some recent progress 
towards understanding the law of the analogue of the 'shape' in this setting has been made in [10, 18]. 

We consider the connection between the tropical RSK correspondence and GL (iV, R)-Whittaker functions 
which we establish in this paper to be quite remarkable and a very interesting topic for further research. 

The outline of the paper is as follows. In Section 2 we introduce the tropical RSK correspondence via 
two equivalent approaches: row insertion procedure and non-intersecting lattice paths. Section 3 provides 
the main set of results of this paper. We state the main algebraic content of the paper in the form of the 
intertwining relation of Proposition 3.4; we define Whittaker functions; we state the paper's main results - 
Theorems 3.6 and 3.8). Within that section we also record the invariant distribution of dynamical tropical 
RSK correspondence, and explain the connection Pitman's 2M — X theorem. Section 4 details how in certain 
scaling limits of our work, one recovers previously discovered results. Proofs of our main results are contained 
in Section 5. 
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Institute. NO'C is partially supported by EPSRC grant EP/I014829/1. TS is partially supported by Na- 
tional Science Foundation grant DMS-1003651 and by the Wisconsin Alumni Research Foundation. NZ is 
partially supported by IRG-246809. 

2. Tropical RSK correspondence 

In this section we describe an extension of Kirillov's 'tropical RSK correspondence' [34] to rectangular 
matrices. We follow mainly the development in [41] but with a slightly different convention for indices. We 
describe first a tropical row insertion procedure, and then expand this into a procedure for inserting a word 
into a triangular array. Repeated insertions create a temporal evolution of the array. In addition to insertion 
into an already existing array we consider insertion into an initially empty array. This latter version will 
have an equivalent description in terms of weights of configurations of lattice paths. 

2.1. Tropical RSK via row insertion. 

Definition 2.1. Let 1 < ^ < A^. Consider two words ^ = {£,g, . . . , ^at) and b = {bi, . . . , &Ar) with strictly 
positive real entries. Geometric row insertion of the word b into the word ^ transforms (^, b) into a new pair 
(^', b') where ^' = (^^, . . . , and b' ~ • ■ • , b'j^). The transformation is notated and defined as follows: 

Ck-^bkiCk-i+^k), e+l<k<N 



b 

(2.1) £, C where 

b' 



b'j, = bkP^, £+l<k<N. 



If ^ = output b' is empty and we write b' = 0. In addition to ^ G (0,oo)^ we admit the case 
^ ~ (1,0,..., 0). This will correspond to row insertion into an initially empty word. With e^'^'' = (1, 0, . . . , 0) 
denoting the first unit fc- vector, the notation and definition are now 

(2.2) AN-i+i) 1^ ^, where i < k < N. 

This is consistent with (2.1) except that output b' is not defined and hence not displayed in the diagram 
above. 



Zl z'l 

0-2 

Z2 ^ Z2 

as 

Z3 z'^ 



Figure 1. Illustration of z' = z <— oi when = 3. Geometric row insertion of the word 
fli = (11170217031) into the triangular array z is defined recursively by insertion of into 
Zi with outputs z[ and Oi+i. After step 3 the process has been exhausted: 03 = (033) has 
one entry and 04 is an empty vector. 



The next step is tropical row insertion of a word into a triangular array. For iV > 1 let Tjv denote the 
set of triangular arrays {zke : 1 < £ < k < N) with entries Zki € (0,oo). The bottom picture of Figure 4 
illustrates an element of T5. (z^i) consists of rows indexed by k and southeast-pointing diagonals indexed 
hyi. 

Definition 2.2. Given z gTn and a word b G (0,oo)^. Geometric row insertion of b into z outputs a new 
triangular array z' G Tn- This procedure is denoted by 

(2.3) z' = z^b 

and it consists of N iterations of the basic row insertion. For 1 < ^ < form words zi = (z«, . . . , zni)- 
Begin by setting ai = b. Then for £ = I, . . . , N recursively apply the map 

ae 

(2.4) ze z; 

a-e+i 

from Definition 2.1, where ag+i ~ a'f. The last output aN+i is empty. The new array z' = (z^^ : 1 < ^ < 
k < N) is formed from the words z^ = (z^^, . . . , z^^). Along the way the procedure constructs an auxiliary 
triangular array a = {aki : 1 < £ < k < N) with diagonals ag = {au, . . . ,aNi)- 

Definition 2.2 oi z' — z b can be summarized by these equations: 

= bk for 1 < fc < iV 

ak+i,i+i = ak+i,e—, tor 1 < £ < A; < A* 

(2.5) ^k+i,e^h.e 

4,f = a.k,e{zk,e + 4-i,f) for 1 < ^ < fc < iV 

4,fc = ak,kZk,k for 1 < fc < A^. 

This procedure is illustrated in Figure 1 when A^ = 3. 

Iteration of the insertion procedure defines a temporal evolution z(n), 71 = 0, 1, 2, . . . , of an array z(n) e 
Tn- This evolution is driven by a semi-infinite matrix d = {dnj '■ n> \,1 < j < N) of positive weights dnj- 
We write c?t^'"l = (dy :l<?<n, l<j< A^) for the matrix of the first n rows of d, and d'"! = {dni, . . . , dniv) 
for the n**^ row of d. The temporal evolution is then defined by successive insertions of rows of d into the 
initial array: given z(0) G Tjv and d, then iteratively for > 1, 

(2.6) z(n) = \z{n - 1) ^ d^''A = [z(0) ^ d^^^ ^ d^l i ^ rfl" 

Figure 2 illustrates. 

Finally we consider the insertion process with an empty initial array. N is still the fixed size parameter 
of the array. Initially z(0) is empty which wc denote by z(0) = 0. The array grows by adding one new 
diagonal Zf at each time. At time n € {1, . . . , A^}, the already existing diagonals zi, . . . , z„_i are updated 
by inserting ai(n) and iterating as in (2.4), and a new diagonal z„ is filled by inserting a„(n) according 
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zjv(3) 





Figure 2. Evolution of the array z{n) in state space Tjv over time n = 0,1,2,... The 
initial state z{0) is on the left edge. The inputs eome from the d-array: ai{n) ~ c?!"! = 

ai(3) 

ziiS) ■■■ 

02(3) 

^2(3) ••• 

03(3) 

^ ^3(3) ••• 

Figure 3. Evolution of the array z(n) started from the empty array z(0) = 0. e^*^' rep- 
resents the word (1,0,..., 0) of length k. By Proposition 2.4 z(n) is equal to the image 
fn,Ar(rf[^'"l) of the weight matrix d)--^'"^ — {ai(i)}"=i under the tropical RSK correspondence. 

to (2.2). Consequently, at time 1 < n < A^, the currently defined array with strictly positive entries is 
z{n) = {zke{n) : l<fc<A^, l<^<fcA n}. We consider the entries {zki{n) : n < i < k < N} undefined. 
At time n = N the array is full, and after time N its evolution continues according to (2.6). The evolution 
of z{n) from z(0) = is illustrated by Figure 3. 

Remark 2.3. Instead of having truncated arrays in the evolution {z{n) : Q < n < from z(0) = 0, we 
could also choose to fill the undefined portion of the array with certain conventions that are consistent with 
the update rules. This would include use of "singular values" and 00. For example, at time < n < A^, 
diagonal Zn+i{n) would equal (1,0,..., 0), in accordance with (2.2). State space Tn would be replaced with 
a larger space that contains and these other partially singular arrays. In this paper we will not use 
these conventions. 

2.2. Tropical RSK via non-intersecting lattice paths. We turn to an alternative definition of the 
evolution in Figure 3 in terms of configurations of non-intersecting lattice paths. As before > 1 is fixed 
and the input of the process is the semi- infinite matrix d = [dij : i > 1, 1 < j < N) of positive real weights. 
For each n > 1 form the n x matrix rfl^'"! = (dy : 1 < « < n, 1 < j < A^). For 1 < ^ < fc < A^ let j. 
denote the set of ^-tuples tt = (tti, . . . , tt^) of non-intersecting lattice paths in 1? such that, for 1 < r < £, tt^ 
is a lattice path from (1, r) to (n, k + r — t). A 'lattice path' only takes unit steps in the coordinate directions 
between nearest-neighbor lattice points of 1? (i.e., up or right); non- intersecting means that paths do not 
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Figure 4. Illustration of the path construction for n = 4 and = 5. Note that the ma- 
trix d is presented in Cartesian coordinates. On the top left is an element tt of H^tv' 
that is an up-right path tt from (1,1) to {n,,N) = (4,5). The weight of this path is 
diid2i<i22<^23'i33(^43C^44(^45- On the top right is an element tt = (711,712) of ^, a pair 
of nonintersecting up-right paths: tti from (1, 1) to {n,N — 1) = (4,4) and from (1,2) 
to (n,7V) = (4,5). The weight of 7r is (^11^21^22^23^33^43^^44) ■ (di2rfi3'^i4rf24rf34d35d45) . 
At the bottom is the array z{n) — {zk,e{n)}i<:(<:k<N , at time n ~ A. We write this as 
z{n) = Pn.N{d). The entry ZN,i{n) = 25^1(4) is equal to the sum of the weights of all paths 
in ^. The product ZM,i{n)zN^2{n) = 25_i(4)z5^2(4) is equal to the sum of the weights of 
all elements in 11^ jy. The rest of the array is determined similarly via (2.9). We regard the 
boldface element Z55(4) as not yet defined at time n = 4 as explained after (2.9). 

touch. The weigiit of an £-tuple 7r = (tti, . . . , tt^) of such paths is 

(2.7) wt{^)=x{ n ^'j- 

'"=1 (jj')eirr 

For 1 < ^ < fc < TV let 

(2.8) Tfc,f(n)= ^ wt{^). 

For 0<n<^<fc<A^ the set of paths ^ is empty and we take the empty sum to equal zero. At i — k 
there is a unique ^-tuple, and in fact we have the equation 

Tfc,f(n) = 5fc/Tfc^„(n) for Q<n<i<k<N 

where 5k,i is the Kronecker delta. For 1 = the right convention turns out to be Tkfl{n) = 1 for 1 < fc < A'^. 
The array z{n) = {zk.i{n) : 1 < k < N, 1 < £ < k A n} is now defined by 

(2.9) Zk,i{n) ■ ■ ■ ZkA'^) = Tk,ein). 

The elements (^Zke{n) : n < £ < k < iV) we regard as undefined, even though strictly speaking one more 
element, namely z„+i.„+i(n), could be consistently defined as 1. In the spirit of Remark 2.3 we could also 
replace the undefined array elements with particular singular values. See Figure 4 for an illustration. 
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We express the mapping (2.9) that defines z{n) from dt^^"! as 

(2.10) z{n) = P„,Ar(d[i^"l). 

We come to the important point from Section 2.2 of [41] that row insertion into an empty array and this 
path construction define the same array z{n). We postpone the proof of this proposition to the end of the 
section. 

Proposition 2.4. Let n,N >1. Set z{n) = P„,jv(d'^'"l) and 

(2.11) z(n) = ^ rfW ^ ^[2] < ^rfM. 

Then z{n) = z{n). 

Let us discuss similarities with the classical RSK correspondence. Array z = P„.Ar(rf'^'"') is the analogue 
of the P-tableau in the usual RSK correspondence. The analogue of the shape of this tableau is the bottom 
vector (zat^i, . . . , ZN,Nf\n)- The analogue of the Q-tableau is the array w = Qn,N{d^^'"'^) = P/v,n((c'^^'"')^), 
where superscript T denotes transpose. It is not difficult to see that the pair (z, w) have the same 'shape', 
that is, {zn,i, ■ • ■ , ZN,N/\n) = {wn,i, • ■ • , Wn,N/\n)- Givcu this constraint, the pair (z, w) can be identified with 
an n X A'' matrix d defined by 



dij — 



Zi+j-i^t, l<i<nAN, l<j<N-i + l 

■WN+n-i-]+i,N-j+i, i<j<N, N-j + l<i<n. 



With this identification, the mapping ^ {z,w) is a bijection from the set ofnx N matrices with strictly 
positive real entries onto itself, which one can refer to as the tropical RSK correspondence. 

Proposition 2.4 shows that, as with the usual RSK correspondence, the 'P-tableau' z{n) = Pn.Af(d'^'"') 
can be defined recursively by inserting the rows of the matrix d one after another. Then w = Qn.N{d) 
plays the role of a 'recording' tableau. With z{k) = Pk,N{d^^'^^) for 1 < fc < 7i, w = Qn,N{d) is given by 
Wk,- = ZN^.{k) for 1 < k < N: this is immediate from the definition. 

Let us introduce some further conventions for the sequel. As before zg = . . . , z^/) is the diagonal 
of z, counting from right to left. The fc"^ row of the array is denoted by z^'^^ = (z^i, . . . , Zkk): and an array 
restricted to a range of rows is denoted by zl^^^l = (zk/ : a < k < b, 1 < i < k) . To discuss evolution of the 
array row by row it is convenient to have notation for spaces of rows. For 1 < k < N the A;*'' row of an array 
in Tn lies in the space of vectors y — {ye ■ 1 < £ < k) with entries in (0, oo). 

As the last item of this section we sketch the proof of Proposition 2.4 from [41]. 

Proof of Proposition 2.4- The connection between z(n) and z(n) goes via the variables Tke{n) and a matrix 
formalism developed in [41]. 

For an Af- vector x = {xi, . . . , xn) define an upper triangular N x N matrix 

H{x) = ^ x^Xi+i ■ ■ ■ XjEij 

l<i<j<N 

where Eij is the N x N matrix with a unique 1 in the (i, j)-position and zeroes elsewhere. For a fixed n, 
define the product 

(2.12) H = i/(dW)i/(fi[2]) . . . iJ(dW). 

A key fact [41, Prop. 1.3] is that the Tki{nys give certain minor determinants of H: Tki{n) = det 
where the superscript specifies the range of rows and the subscript the range of columns in the minor. 

On the other hand, the row insertion procedure can be encoded with i/-type matrices. Let 1 < m < A 
and introduce this further definition for an (A^ — m + l)-vector x = (a;,„, Xm+i, ■ ■ ■ ,xn)'- 

(2.13) Hr,i{x) ^ ^ Pi,i + ^ x^Xi+i ■ ■ ■ XjEij. 

l<2<m m<i<.j<N 

In particular, Hi{x) = H{x). 

Set AI = n A N. With H as in (2.12) above, consider the equation 

(2.14) H = HMimi) ■ ■ ■ H2{m)Hi{m) 



for unknown vectors rji = {r]k/)e<k<N , i = 1,2,..., M. This equation is uniquely solved by [41, Thru. 2.4] 

(2.15) Tji i = — , r^k,i = 1—\ TT for £ < fc < TV. 

Equation (2.14) encodes the row insertion procedure but in different variables [41, eqn. (2.38)-(2.40)]. 
Namely, the 77-variables are the ratios of the z-variablcs defined by (2.11): 

(2.16) for 1 < £ < M: = ^^/(n) , = iov £ < k < N . 

Zk-i,i[n) 

Combining (2.15)-(2.16) for z with (2.9) for z gives 

(2.17) 5fc,,(n)=77„---?7fc,, = ^Mi^=2 („) foTl<e<M,i<k<N. □ 



3. Tropical RSK with random input 

Given an initial (possibly random) state z(0) € Tat and a weight matrix d composed of independent 
random rows with z(0) and d independent. Proposition 2.4 shows that z{n) — Pn,jv('^'^'"') has the 
structure of a Markov process with time parameter n. The exact form of the transition kernel depends on 
the distribution of the rows dl'l and can be explicitly written down by appealing to the recursion of Definition 
2.1. We do this for a particular solvable distribution on the elements dij that we now introduce. 

Definition 3.1. Let be a positive real. A random variable X has inverse- gamma distribution with param- 
eter > if it is supported on the positive reals where it has distribution 

(3.1) P(X e da;) = ^2;-''-iexp|-i|da;. 
We abbreviate this X - 1(6*). 

Definition 3.2. An inverse-gamma weight matrix, with respect to a parameter matrix 7 = {"^ij > : i > 
1,1 < i < N), is a matrix of positive weights {dij : i > 1, 1 < .7 < N) such that the entries are independent 
random variables and dij ^ r^^(7,;.j). We call a parameter matrix 7 solvable if 7,;.^ = 0, + Oj > for real 
parameters {Oi : i > 1) and {6j : 1 < j < N). In this case we also refer to the associated weight matrix 
as solvable. Column n of the parameter matrix 7 is denoted by 71"] = {'jn,j)i<j<N ■ We denote the vector 
9 ^ {9j :!< j < N) for later use. 

The intertwining properties we discuss in Section 3.1 are the reason we restrict to inverse-gamma dis- 
tributed weights and the reason for the particular form 7^.^ = 0i + 6j. 

The transition kernel for the Markov chain z{n) on the state space T^v resulting from applying the tropical 
RSK correspondence to a solvable inverse-gamma weight matrix is denoted by n;^„] (z, di). This represents 
the time n transition z{n — 1) — > z{n). 

To explicitly state this kernel, it turns out useful to exploit another structural property of the image of 
the tropical RSK with independent weights dn^k- the rows of the array z{n) form a Markov chain (indexed 
from top to bottom) with respect to adding columns to the weight matrix d. For this purpose let us denote 
z = z{n — 1) and z = z{n). We begin at the top of the array. The singleton top row (denoted by y) of z is 
updated at time n by the transition kernel 

(3.2) p;H(y,#) = r(7„.0-(^)"%xp|-^j?l. 

This simply encodes y ~ dn.iy with dn,i ^ r^^(7„^i). 

Now we move down along the rows of the array (recall z^'^^^l = {zk-i/)i<i<k-i and likewise for z^''^^^). 
Given both the initial and updated row fc - 1, (z['=-il,z['=-il), from (2.5) we read off the time n rule for 
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updating row k from = {zk^e)i<i<k to zl'"' = {zkd)i<e<k- 

The new input weight is ak,i ~ dn.k ^ r~^(7„.fe), and the equations are 

' Zk,i = ak,iizk,i + Zk-is) 



(3.3) 



Zk.e-iZk-i.e-i Zkj + Zk-ij r n ^ /, ^ , i 
Zk,e = — '■ • '■ r '■ for 2 < i: < — 1, 

Zk-i,e-i Zk,e-i + Zk-i,e-i 

Zk,kZk,k-lZk-l,k-l 

Zk,k ~ 



{Zk,k-1 + Zk-l,k-l)zk-l,k-l 



We encode this step in a kernel L^^ ^((•^''^ ^1, zl*-'!; z^*^ -^l),dz['^l) from Yfe_i x x Yk-i into Y^, given as 
follows in terms of the integral of a bounded Borel test function g on Y^ : 



(3.4) 



J {0,00)" ' J (0,00) yi \ VI J I yi J 

yi-ixi-i ye + Xi \ ykVk-iXk-i 



X 9\ Vi 



xe-i ye-i + ^£-1 J 2<i<k~i ' ^k-iivk-i + Xk-i) 



Now we can write down the kernel for the evolution of the array. The kernel n;^^„] (z, dz) for the transition 
from z — z{n — 1) to z = z(n) on the space T^r is defined inductively on N . For = 1 set n^[„] = P^[„j, 
and for > 2 

(3.5) n'^,„,{z^''''\d~z^''''^)=n^,-,\z^^'''~'\d~z^''''-'^)L^^^^^ 

3.1. Intertwining relation. Let us first recall a well-known criterion for a function of a Markov chain to 
retain the Markov property (see, for example, [ol]). Consider a measurable transformation cf) : T ^ S where 
(T, T) and {S, S) are measurable spaces. Given Markov transition kernels n„ on T one forms a Markov 
process z(n) which has a given initial distribution z(0) and transitions between z{n — 1) to z(n) via n„. The 
process {z(n)}„>o is Markovian with respect to its own filtration. The question is, under what conditions is 
y{n) = (j){z(n)) Markovian under its own filtration a{y{0), . . . , y(n)}, and what is the associated transition 
kernel Pn'-S ^ S7 

In order to answer this question we introduce an intertwining kernel K : S ^ T. (The reason P and 
K have bars is that one often initially deals with unnormalized kernels and then normalizes them to be 
probability measures in their second variable.) 

Proposition 3.3. Assume that there exist Pn and K (which does not depend on n) satisfying 

(i) for all yes, K{y,cj>-^{y)) ^ I, 
ill) for all n, KU„ = P„K . 

Then, for any initial (possibly random) state y'^ e S, if one initializes the Markov chain z(ri) with z(0) 
distributed according to the measure K(y'^,-), one has these properties: 

(i) For all y e S and all bounded Borel functions f on T , 

E[/(z(n))|y(0), ...,y{n- l),y{n) = y] = {Kf){y). 

(ii) The process y{n) = (l){z{n)) is Markov in its own filtration a{y{0), . . . , y{n)} with transition kernel 

Pn- 

Return to tropical RSK with the solvable inverse- gamma weight matrix. Take T above to be T^r and S to 
be Y^r. Let </> : T^r Yjv project z G T^f onto its bottom row. It is easier to first introduce unnormalized 
kernels and prove intertwining, and then to normalize them to apply the above results. 

Define a time n positive kernel on Y^r by 

no«p{-^}n(ru...-(|)'"=p{-|}f)- 
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Define a positive intertwining kernel from Yjv to T^r by 

(3.7) K-iy,dz)= n f^)'"^'^xpf-^-f^±M±iWfj<5^^(rf,,,) 

where 5y{dzij) is the Dirac delta measure at y. Kg is the identity kernel. Observe that Kg only depends 
on the column parameters 9 which do not change with the time index n. Definition 3.2 stipulated the form 
7ij" — di + 6j for the solvable parameter matrix in order to make the intertwining work. The intertwining 
itself does also work with a time-dependent if-kernel. But in our case computations reveal that application 
of Proposition 3.3 requires a time-independent K and we are not permitted any more general (ji.j). 

P^[„] is a stochastic kernel (i.e., normalized to have measure one) and represents the update zi.i = ai.izi.i. 
But for N > 2, P^[„] is substochastic. This is evident because the second product is the stochastic kernel of 
independent inverse-gamma distributed multiplicative jumps, while the first product is a killing potential. A 
Doob /i-transform (or ground-state transform) will suffice to renormalize this kernel as well as the intertwining 
kernel. 

The main algebraic content of the integrability or solvability of tropical RSK is: 
Proposition 3.4. The following intertwining relation holds at all times n > 1: 

(3.8) P^^iAT = <n^.i, 
where both sides are operators from Yn to Tn- 

This is proved in Section 5.1. 

The kernels above arc not normalized. However, using the intertwining relation it is now simple to 
determine the necessary normalizing functions. For y g Y^v, define 

(3.9) w^{y)^ I K^iy,dz). 
Integrating the intertwining (3.8) yields the eigenfunction relation 



(3.10) P^.^.uj^ = 



jN N _ N 

Thus we can define a stochastic kernel on Yat by 

(3.11) Pt,(y,dy) = ^P^^,(y,dy) 
and from Y^v to Tat by 

(3.12) K^{y, dz) = -^K^{y, dz). 

wg (y) 

The kernel K^{y,dz) should be interpreted as the distribution of the pattern in T^y conditioned on the 
bottom row z^^l being equal to y. From the previous proposition follows: 

Corollary 3.5. The following intertwining relation holds at all times n>l: 

(3.13) P^l,K^=K^n^,„„ 
where both sides are operators from to T^. 

3.2. Whittaker functions. For the next stage, note that the kernels above remain perfectly well-defined 
if we allow parameter vector 9 to be complex. The probabilistic meanings arc lost but the intertwining 
continues to work. For y G (0, oo)^ and A e C^, define 

N 

(3.14) (y, dz) = n (y, dz) 
and 

(3.15) ^^(y)-^ f M^{y,dz) = f[yr^^w^{y). 
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The functions well-defined for any A € C^, are class-one GL(iV, M)-Whittaker functions (in multi- 
plicative variables). They arise in various contexts: they are eigenfunctions of the quantum Toda lattice 
(when expressed in additive variables Xi = log yi) and can be represented as particular matrix elements of 
infinite-dimensional representations of 3l{N) [3()]; they also arise in the harmonic analysis of automorphic 
forms on GL{N,M.) [14]. The integral representation (3.15) is due to Givental [28]. It is known [53, 33] that 
the integral transform 

(3.16) /(A)=/ n— 

defines an isometry of L2((0, od)^, J|- dyi/yi) onto L2^™(iM^, sjv(A)(iA), where Lj^™ is the space of L'^ 
functions which are symmetric in their variables, l = \/— 1 is the imaginary unit and 

(3-17) ^^^W-7^^^nr(A,-A.)-\ 
The inversion formula is 

(3.18) 5(2/) = / giXMiy)sNi\)dX. 
In particular, the Plancherel formula 

(3.19) / f(y)J(f)Y[^= f f(X)J{X)s^{X)dX 

holds for functions f,g & L^((0, oo)^,n dyi/yi). 

We also have the Whittaker integral identity [15, 58, 27], for s > and A, e C^, 

(3.20) / e-^y^^^{y)^:^{y) {[^ ^ [] ^(-A. - 
Using ^'e(y) = '^^g{y'), where y[ = yj;\i^i, this is equivalent to 

(3.21) / e-^y-^'^^iy)^^iy) = J| r(A. + .,). 
Note that if, for z G Tjv, we define Xi{z) by 

k k 
Y[x,{z)=Y[zk,^, k=l,...,N, 

i=l i=l 

then 

N 

(3.22) Mg^{y,dz) = l[x,{z)-'^ {y,dz), 

i=l 

where 

(3.23) M^{y,dz)= ^ e.p ( - ^ - ^ f[Sy,{dz^,). 

\ Zk+l f Zh p / ZhP -^^ 

l<e<k<N ^ t+i.t / n,i 

For n > and i — 1, . . . , N write Xi{n) — Xi{z{n)). Then Xi{n) is a multiplicative random walk: 

(3.24) Xi{n) = [ Y[ 

(0). 

\m=l I 
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3.3. Main theorems. We are prepared to state the two main theorems of the paper. They are proved in 
Sections 5.2 and 5.3, respectively. The first result is concerned with the solvability of the </> projection of the 
z{n) Markov chain corresponding to the recursive system in equation (2.5). 

Theorem 3.6. Fix a solvable inverse-gamma weight matrix defined in terms of parameters {9m '■ m > I) 
and {9j : 1 < i < N) and assume (without loss of generality) that 9j < < 9m for all j,rn. Let y{Q) be a 
random or deterministic initial state in Yjv and let the initial distribution of z{0) be K^{y{0), ■ ). 

(i) The sequence of random variables y{n) = (j)(z{n))^n >0, is a Markov chain with respect to its own 
filtration, with state space Yjv, initial state y{0) and time n transition kernel P^[„]- 

(ii) For a bounded Borel function f on Tn and y € Y^r 



(3.25) E[fizin)) \ 2/(0), ...,y{n- l),y{n) = y] = j K^{y, dz) f{z) 
(iii) For A e 

(3.26) E [ [] .T,(n)-^- 1 y(0), . . . , y(n - l),y{n) 



N ,T,JV 

y 

1=1 

(iv) For an initial state G Y^r and time > 1, let /i^ {y^,dy) denote the probability distribution of the 
time-n state y{n). Then for all X E lM.^ , 

„»,..0 J.,> _ *A (»°) A A r(«.n + A.) 



(3.28) 



Moreover, for any continuous, compactly supported function f on (0, oo)"'^ we have 



I iM{y\dy)=f dxs,{x)'§^(f /(y)<(y)*^!.(2/)n-) nnj^l^- 

Remark 3.7. Note that, given part (i) of the above theorem, (3.27) shows that one may diagonalize the 
transition kernel via the eigenfunction equation 



(3-29) ^.H^- n 



^^_/^TT m. + A,-) \^f 



■ j= 

which can also be seen directly from the intertwining relation (3.8), cf. (3.10). By applying the completeness 
relation resulting from the isomctry, this identity characterizes the transition kernel. 

Next we specialize the above result to the Markov chain y{n) that comes from the evolving shape of the 
tropical RSK array P„,Ar(d[^'"l) of (2.10). This is the case of the empty initial array. We can capture this 
situation by taking a somewhat delicate limit of the initial state y^. This is our second main result. 

Theorem 3.8. Fix a solvable inverse-gamma weight matrix defined in terms of parameters {9m ■ m > 1) 
and {9j : 1 < J < N) and assume (without loss of generality) that 9j < < 9m for all j, m. 
Consider the array p = {pk.e)i<e<k<N with 

fk-l k-l 

(3.30) pk = {pk,i)i<i<k ^ , — 1,. 

for l<k<N. Let y^^^' ^ {^~^'^'''')i<i<n "''^'^ n>N. Then 

(i) As M — > oo, the probability distribution p^{y^'^\dy) converges weakly to a distribution p,^{dy) 
characterized by 

/(2/)M^(rfy)= / dXsN{X) ' 



/ /(.)*.(.)*-.(.)nv nn 

\'^^+ i / m=l 1=1 



r{9m + AQ 
r{9, + 9m) ' 



for any continuous, compactly supported function f on (0, oo) 
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(ii) The Laplace transform of the projection of /i^ (dy) on the first coordinate is given by 

N 



(3.31) / e-^y^f,^idy)= f dX s^^^^'^^'^^ [] r(A. ^ H 11 



SAr(A), 



r(0i + 9m) 

where the poles of the functions T{Xi — 9j) and T{\i + Om) are not encountered due to the assumed 
condition that d„i > for all m and 9j < for all j . 

(iii) The measure ji^ (dy) is the distribution of the bottom row y{n), given that the process begins with 
the empty array. In particular, the distribution of the partition function z^r i(ri) = X^Tren^ wt('K) 

is the marginal distribution of fi^ (dy) on the first coordinate yi, and hence uniquely characterized 
by (3.31). 

(iv) The distribution of Pn,N{d^^'^^) is the measure in dz given by 

fi:!idy)K^{y,dz). 

(v) When n ~ N we have the following simplification: 

N N N , 

t^%{dy)=l[ l[r{0, + Om)-'e~y-'^e{y)^Ml[—- 

m=l i=l i=l 

In particular, for s > 0, 

e-^^«Vi^(rfy) = (i + .'^)-^'-^''+^'\ 

(0,oo)" 

that is, the random variable zn,n{N) is inverse gamma distributed with parameter X]i!=i(^« + ^«)- 
The distribution of zn^n{N) can also be seen from (2.9). 

Observe that the condition of n > is not restrictive when it comes to computing the Laplace transform 
in part (ii) of the above theorem. Indeed, if one wishes to compute the Laplace transform for n < N then it 
suffices to transpose the parameter matrix and switch the role of n and N. The distribution of the coordinate 
yi is unchanged by this procedure, and now the above corollary applies. 



3.4. Pitman's 2M — X theorem. Theorem 3.6 can be regarded as a variant of Pitman's '2M — X theorem', 
which states that, if Xt is a standard one-dimensional Brownian motion and Mt = maxs<t X^, then 2Mt — Xt 
is a three-dimensional Bessel process. This theorem has vast generalizations [G, 8, 9, 13, 23, 35, 39, 42, 43, 44, 
45, 48], many of which have been obtained via various analogues the 'Burke property' discussed in Remark 
3.10 below. All of these can be regarded as variations of the statement that the stochastic evolution of the 
shape of the tableaux, obtained when applying variants of the RSK algorithm to random input data, has the 
Markov property. The first 'geometric' or 'positive temperature' version of this statement was discovered 
by Matsumoto and Yor [39], who showed that, for Xt as above, the process log/oe2'^=-^*ds, t > IS a 
diffusion on M with infinitesimal generator given by 

1 ^ / d ^^^j^ (e~^)^ ^ 

2 dx"^ \dx J dx' 

where Kq is the Macdonald function (with index 0). A multi-dimensional version of this theorem of Mat- 
sumoto and Yor is given in [42] , which can be regarded as a particular specialization (scaling limit) of the 
main result in the present paper. It is also proved via an intertwining relation and is closely related to 
the quantum Toda lattice. The corresponding directed polymer model is defined on the semi-lattice Z x R. 
Both models feature the G'L(A^, M)-Whittaker functions in an essential way; the relation between them is 
analogous to the relation between the Gaussian and Laguerre unitary ensembles in random matrix theory. 
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3.5. Invariant distributions. The Markov process defined by row insertion with a solvable inverse gamma 
parameter matrix turns out to have nice invariant distributions. The z-array itself cannot have an invariant 
distribution: for example, Zi^i(n) = dn.i ■ ■ ■ (ii,izi,i(0) evolves as a multiplicative random walk. Instead, we 
look at ratios of z- values. 

Fix > 1. For an array z G Tn define the array i] = {rike)i<e<k<N of ratios by 

Vk.t = , 1 < ^ < A: < iV. 

Zk-i,e 

The Markov process z{n) then defines another random process r]{n) = {r]k£{n))i<i^k<N by rjk.iin) = 
Zk,i(ji) / Zk-i,i(n). Denote again diagonals by rii(n) = {'r]ki{n))i^k<N for 1 < ^ < A^. This new process 
r](n) will also be a Markov chain. 

Theorem 3.9. Let z{n) evolve on the space T^r according to the Markovian dynamics governed by a solvable 
inverse-gamma weight matrix with parameters = 9i + 9j, as specified by the transition kernels in (3.5). 

(a) The process rj{n) is a Markov chain in its own filtration. 

(b) Let 1 < j < N . Assume 9i < 62 < ■ ■ ■ < Oj < min{6'j-f 1, . . . , 6^}. Then the process (771(7?,), . . . , r/j{n)) 
has an invariant distribution where the variables {rj^e ■l<i<j,i<k< N} are independent with marginal 
distributions rjki ~ T~^{9k — 9i). If the process is started with this distribution, then the following statement 
holds for all times n> \: the variables {r]ke{n) '■l<£<j,(<k< N} U {zNi{m) / Zffe(m — 1) : 1 < m < 
71, 1 < £ < j} are independent with marginals r]ke{n) ^ T^^ {9k — 9^) and ZNi{m) / zjqi{m—l) ^ T^^{9m + 9i). 

Remark 3.10. Theorem 3.9 is an extension of a result of [0(1] for the directed polymer with inverse-gamma 
weights. This could be called a 'Burke property' by analogy with the Burke theorem (also known as the 
'output theorem') of M/M/1 queues. According to the Burke theorem, for a reversible queue the number 
of customers in the system at time t is independent of the departure process up to time t. This notion has 
an analogy in models with random weight matrices, and it has been used in the past to derive exact limit 
shapes [■j4, 'm] and fluctuation exponents [4, 5, IG]. In fact, it was this property that led us to investigate 
the solvability of tropical RSK for inverse-gamma weight matrices. The analogous Burke property was found 
earlier in the Brownian polymer model [47] and was used to derive fluctuation exponents for that model in 
[57]. 

Theorem 3.9 is proved via (2.4) that represents a transition of the Markov process z{n) in terms of a 
sequence of tropical row insertions. For this purpose we reformulate the row insertion step in terms of the 
ratios. In Definition 2.1 with fixed 1 < ^ < the inputs of the row insertion were ^ = {^i, . . . and 
b = (6f , . . . , ^Af), and the outputs ^' = (^^, . . . , ^^) and b' = (fef_|_i, . . . , &^). Define now r^k = £,k/£,k-i and 
v'k = Ck/Ck-i for i < k < N, and also auxiliary variables (k = Ck/^k ioi i < k < N. The words are 
V = ive+i, ■■■,Vn) and 77' = (77^^^, . . . , ry^). 

Lemma 3.11. Fix integers 1 < £ < N . In terms of the new variables, tropical row insertion transforms 
{rj, b) into (77', b') via the following equations. Set first Q = be, and then inductively for k — £ + 1, . . . , N : 

(3.32) ^^ = 5fc(l + -^), a = bk(l + ^), and bi=(-^ + 

^ Qk-i^ ^ Ilk ' ^ Cfe-i '7fe ' 

Next the row insertion step with random input. 

Lemma 3.12. Fix integers 1 < £ < N . Let ae+i, . . . ,aN, f3e,...,PN be positive reals that satisfy jSk = 
f^e + ctk for £ < k < N . Assume that the random variables {rjk : £ < k < N} U {bk : £ < k < N} 
are independent with marginal distributions rjk ^ V~^(ak) and bk ^ V~^(j3k). Then the random variables 
{^fc'^fc ■ £ < k < N} U {Cn} also independent with marginal distributions rj'/. ^ T^^{ak), b'j^ T~^{l3k), 
and Cn ^ (/?<>). 

Proof. From the assumptions and by definition, Q = bi ^ T^^{j5i) and this variable is independent of 
{77^+1, . . . , 777V, bi+i, . . . , bisi}. Use equations (3.32) to prove, inductively on 771 = £ -f 1, . . . , A^, that random 
variables {rj'g^-^, . . . ,ri'^,b'^j^-^, . . . ,b'^,C,Tn} are independent, independent of {77„i+i,..., 7?Ar, 6„i+i, . . . , ^at}, 
and their marginal distributions are 77^ ^ T~^{ak), 6'^, ~ r~^(/3fe) and ^ T~^{(3i). An induction step 
is achieved by applying (3.32) to the triple (Cm, '7771+17 &m+i) to produce the new triple (Cm+i, ^7™+!: ^in+i)- 
Note that the parameter of Cm does not change with m. The case m = N gives the lemma. □ 
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Proof of Theorem 3.9. (a) That rj{n) is itself a Markov process follows from the fact that from (3.32) we can 
build autonomous equations for this evolution. 

(b) It suffices to show that the last claim is preserved by a step of the evolution. Consider the time n 
transition from state ?/ = ri{n — 1) to state ry' = rj{n). The input weights are ai ~ (an, . . . ,aNi) = c^'"' = 
^n) with dn,k r ^(7n,fe), and this also defines the first diagonal ai of the auxiliary array in 
Def. 2.2. Assume that the variables {r]k£ ■l<i<j,i<k< N}U{zNi{m)/zNi{m-' l):l<m<n — 1, 1< 
^ ^ j} are independent with marginals r]M ^ T~^{0k — 6() and ZNi{m)/ ZNi{m — 1) ~ ^"^{Om + (^e)- Let 
(^Ni = ZNi{n) / ZMi{n — 1) denote ratios defined along the transition process. 

We prove the following statement inductively over ^ = 1, . . . , j. 

The variables {rj'i, . . . ,ri[, a^+i, Cni, ■ ■ ■ i Cni} are independent and independent of 
(3.33) {Vi+i^---iVj}^{zNi{rri)/zNi{rn-l):l<rn<n-l,l<i<j}, and 

they have marginals rj'^i ^ r~^(6'fc - Oi), ak/+i ~ r~^(^„ + Ok) and C,ni ^ r"^(^„ + 6*,;). 

The case (. = j completes the proof. 

In the first row insertion step apply Lemma 3.12 with £ = 1 and inputs rji = (772,1, • • ■ iVn,i) and b = 
(an, . . . , otvi)- Now ak = 0k — ^1 and = jn,k = (^n + ^fc- According to the lemma, the outputs 
Vi = {v'2 It ■ ■ jV'n 1)' b' ^ 02 ~ (a22, • ■ • 7 aAr2) and Ctvi are independent and they have the correct marginal 
distributions: 77^^ - T-^idk - Oi), ak2 r~i(^„ + 0k) and (ni ^ ^'H^^n + Oi). This gives (3.33) for £ = 1. 

In the general step, assuming (3.33) for (—1, Lemma 3.12 is applied to inputs 77^ and ae, with ak = 0k — di 
and j3k = ^n,k = + Ok. The outputs 77^, a^+i, C,n,i have the right properties and the validity of (3.33) is 
exteded to ^. □ 



4. Degenerations to known results 
We detail rescalings of the inverse-gamma polymer which recover known results. 

4.1. Directed last passage percolation and the Laguerre Unitary Ensemble. Fix a solvable param- 
eter matrix 7 — {jij >0:7>l,l<j< N) such that jij = 0i + 0j. Consider a family (indexed by e > 0) of 
solvable inverse-gamma weight matrices d"^ = (df ^ : i > 1, 1 < J < N) such that the entries are independent 
random variables and df ^ ^ T~^{eji,j). Keeping track of the e write z| ^(n) as the elements of z^{n) - the 
image of the weight matrix {df ^ : 1 < « < ti, 1 < j < N) under the tropical RSK correspondence. Write 
F^n) = {Fy{n) ■.l<l<k<N) where F^i,{n) = e logz|_^(7^). 

With respect to the same solvable parameter matrix consider a weight matrix w = {wij : 7 > 1, 1 < j < N) 
such that the entries are independent random variables and Wij ~ Exp(7ij ) (an exponential random variable 
with rate 7^.^, or equivalently mean (7i.j)^^). The classical RSK correspondence maps the weight matrix 
777 to a pair of Young tableaux {P,Q) and is defined analogously but with the (-I-, x) algebra replaced by 
(max, +). We focus on the P-tableaux and writing it in terms of a Gelfand-Zetlin pattern (a triangular 
array with interlacing). Recall the notation for non-intersecting paths 11^ from Section 2. The weight of 
an ^-tuple tt of paths is now 

e 

(4.1) wt{TT)=Y^ E ^'J- 

r=l 

Define an array ^(77) = {Lk.i{n) : 1 < k < N, I < I < k An} hy 

(4.2) Lfc 1 (77.) + • • • + La:,^ (77) = max wtM. 

L{n) is the Gelfand-Zetlin pattern version of the P-tableau of the image of {wij : 1 < 7 < n, 1 < j < A^) 
under the RSK correspondence. As in the case of the tropical RSK correspondence, when 77 < A^ it is 
necessary to leave some entries of L{n) undefined, or populate them with singular values (see Remark 2.3). 

Proposition 4.1. The n-indexed process iF^{n))^_^~^^ converges in law, as e — > 0, to the n-indexed process 
iL{n))n>,- 
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Proof. This hinges on two observations. The first is that as e goes to zero, elogdf j converges in distribution 
to Wij. Hence by independence the whole array {elogdf j} converges in distribution to w and (by the 
continuous mapping theorem) the random vectors 

(4-3) (e E '^^d^A =^ (E E 

The second observation is that on R™, the function fe{x) = £^ogYlT=i ^^'^^ converges uniformly, as £ — ^ 0, 
to the function fo{x) = max(xi : 1 < i < m). The process (-F'^('t-))„>o is formed by applying an array of 
functions of the type to the elements of the vector on the left of (4.3). Combining the uniform convergence 
of functions of this type with the convergence in distribution in (4.3) , the claimed convergence of the processes 
follows. □ 

Let us now briefly recall the connections between bottom row of Lin) and the eigenvalues of the Laguerre 
Unitary Ensemble (LUE). Consider an array (A^.j '■ ^ < i < N ^ j > 1) oi independent complex zero-mean 
Gaussian distributed random variables with variance {"ijji)"^ ■ We have changed the order of i and j since 
we now set A{n) = {Ai^j '■l<i<N,l<j<n) and treat A{n) as a matrix with (row, col) notation. 
Set M{n) — A{n)A{n)* the N x N generalized Wishart random matrix, and define for each n, a vector of 
ordered (largest to smallest) eigenvalues of M{n): \{n) = (Ai(n), . . . , AAr(n)). 

When 7i.j = 1 for all Johansson [-31] showed that for fixed n, \i{n) Lpf i{n). This was strengthened 

by [22] to show that for the same parameters as in [31] and for n > fixed, the vector (Ai(n))^j^ =■* 
{LN^i{n))f^^. In [20, 25] equality in law was shown for the processes (^i("-))„>i and {LN,i{n))^^y^. 

Turning to the general case for the parameters jij, [12] proved that for n fixed, Ai(n) Lm,i{ji) and 
conjectured equality of the corresponding processes in n. This was then proved in [21], whose methods 
(combining Theorem 3.1 and Lemma 4.1) show equality in law of the process (A(n))^^-|^ and (-^Af('T-))„>i 

where L^in) — {L]y,i{n))^^^ and where only non-zero eigenvalucs/non-singular entries of L are considered. 

Combining Proposition 4.1 with the above discussion one sees that the logarithm of the bottom row of the 
image of a solvable weight matrix under the tropical RSK correspondence are analogous to the eigenvalues of 
the LUE ensemble. This connection can be seen from the integral formulas we have derived in (i) of Theorem 
3.8. We perform point-wise asymptotics to demonstrate this connection we have proved above. Performing 
the change of variables to that formula given hy yi i-^ e'^ ^% 6* i—> ed and A i-> eA we find that the measure 
for {F^j{n))^^^ = {xi,..., xn) is given by 

-SAr(eA). 



7=1 -^'R" rn=lT=l 



r(£(0, + e^)) 

We may now evaluate the e — !• asymptotics: For the Gamma functions we employ the expansion near 
zero; for the measure Sn we employ the Euler Gamma reflection formula; for the Whittaker functions we 
can write the Whittaker functions in additive variables (and perform a sign change) ip\{x) = ^^^{e^), and 
then use the fact (sec for instance [33, 42]) to see 



lim £^^f.(e-x) = ^-(x) = ^-,...,,„(.) = (-1)^ det(e>-05^r 



2 

e^O - " - ' - - h(\) 



where h{\) = ni<i<j<Ar('^j ~ ^i) is the Vandermonde determinant. 

The expansion of the Gamma function near zero shows that as e goes to zero. 

Likewise, the Euler Gamma reflection formula and the fact that sin(7ra;)/(7ra;) — )■ 1 as x — > yields 
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where h{X) = ni<i<j<Ar(-^j ~ ^i) ^^e Vandermonde determinant. 

Putting these asymptotics together we have that the Umit of our measure is given by 



N 



1 



1=1 



N 



dAdet(e^-05=iMA)nn- 



h{e) {2ni)^m 
For simpUeity let us assume that N = n. Then this can be rewritten as 

N 

Zn,n 



m=l i=l ""I 



where we have 



D{x,e) 



{2m 



— det(e-«'"Of,=i^(^,^) TT'^^i 
„ / dAdet(e^-05^iM^)MA) IT 11^— V' 



where the integrals are along lines parallel to the imaginary axis and to the right of the poles, and where 



Z 



N,N 



h{e)h{e) 



1 



= det(— 



n 



N 



An applieation of the Residue Theorem provides that 



N 



and henee our measure is simply 

1 



AT 



det(e-'''^05=i det(e-^'^05=i H dx^U>o 



which coincides exactly with the formula given directly for last passage percolation and for the generalized 
Wishart ensemble in [12]. 

4.2. Semi-discrete directed polymer in a Brownian environment. We will now indicate how an 
appropriate scaling of the solvable parameter matrix (7ij ) can be used to recover the semi-discrete directed 
polymer in a Brownian environment studied in [47, 40, 57, 42], as a scaling limit of zjvi(n). In particular, 
we will consider the weight matrix (dij) with solvable parameter matrix jij = n and we will let n tend to 
infinity. We would first need some facts about the digamma and the trigamma functions, which are defined 
as 



(4.4) 



^o{0) ^ -^logm and M/i(0) = ^logr(0). 



In particular, we have that for 9 large 

(4.5) vl/o(0)^log0-l+o(i), vl/i(0)^i + o(i). 

We also notice that for an inverse Gamma random variable d with parameter it holds that 

*o(6') = -E[logd] and 4'i(6') = Var(logd). 
We write logzjvi('T.) in terms of the weights dij (iid inverse Gamma random variables with parameter n) as 

N 



logzjvi(n-) = log Y "^^P 

l<.ii ■ ■ ■<ij^—n 



exp 



logn-(^-i) J2 

l<.ii ■■■<.ij\j—n 

+{N ~ l)logn- {n + N - l)*o("-) 
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j = l ^ ij-l<i<ij 



Using (4.5) we have that n^'i(n) -> 1 and {N — l)logn — {n + N — l)'I'o(n) = —nlogn + | + o(l), as n 
tends to infinity. It is now an easy consequence of Donsker's invariance principie and the Riemann integrai 
definition that 

log (n"z^i(n)) - 1 =^ log / . . . / eE£.(B'(tO-B'(t.-.))di^ . ..dt„_,, 

^ J Jo<ti<-<t„-i<l 

for n tending to infinity, with independent, standard Brownian motions. This is the directed polymer 

model studied in [47, 40, 57, 42]. Along the same lines it follows that the whole pattern {zki{n)) obeys 
similar scaling limits. 



5. Proof of main results 

5.1. Proof of Proposition 3.4. We prove the intertwining relation (3.13) in two steps. We use induction on 
N to show that (3.8) is valid. The inductive step will also reveal an inductive property of the eigcnfunctions 

The case = 1 in (3.8) is immediate since n^[„j — P^[n] and A'g^ is the identity operator. 
A supporting step of the proof is an intertwining that involves only two rows of the array z. For this 
purpose introduce two further kernels for 2 < A: < A^ < n: a time n kernel on Y^._i x Yj. by 

(5.1) {{z^-\z^),dz'^~\dz'^) = P^^-^{z'^-\dS'-^) L';^^ J{z'^-\ 2^ z'^-i), dz'^) 

and a time- homogeneous kernel from to Y^-i by 



e=i 



ye 



k-l / ST k-1 



ye xi 



dxf 



1=1 



Lemma 5.1. At every time n > 1 and for all y G Y^v, we have the following equality of measures on 
Yat-i X Yn, in terms integration variables {dz^~^,dz^): 



(5.3) Pj^Uy,dz'')A'^{z'\dz''-')= / A'^{y,dx)R'^,^,{ix,y),dz^-\dz''). 

With (5.2) we can give this alternative representation to the intertwining kernel (3.7) from Yjv to Tjv: 

(5.4) K{y,dz) = Syidz"") n A^-^-(z^-^dz^-^-i). 

Postpone the proof of Lemma 5.1 for a moment. With (5.4) and Lemma 5.1 we can complete the proof of 
Proposition 3.4 by checking the induction step. Assume (3.8) for A^ — 1. 

TV — 2 

/ P^Uy,dy)K^{y,dz'-^) = P^Uv^dz^) [] A^~' (z'^-^ ,dz^~^'') 
by Lemma 5.1 

„ N-2 

= / A^{y,dx)R^,„,{{x,y),dz''-\dz'') H A^-\z''-\dz''~^'') 



by definition of R^[„] 



I A^{y,dx)P'i-,\x,dz^-') n A,^-^"(z^-^dz 

xL^^J{x,y;z^-'),dz^) 
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by definition of Kn and tlie induction assumption 

A^{y,dx) 



by noting that z^^^ — x under K^^~^ ^{x, dz^-'^~^), and by definition of 11 



Ae{y.dx) 



K. 



N-l 



Kl^iy,dz^-'')U':^,„,{z^'^\dz^-^'). 
This checks (3.8) for N. 

Proof of Lemma 5.1. Take a test function g on Y^v-i x Yjv, and collect and rearrange all the factors in the 
integral against the kernel on the right-hand side of (5.3): 



Kiy,dx) 



R!^W]{{i,y),dx,dy)g{x,y) 



TTr(7«.j) r / dx dyi dx 



X exp 



N-l 



N-2 - , _ 



e=i 



X£ ^ ye+i _|_ X£ 
, xe xi 

N-l 

x(yi+ii)^-"yr-""' • n 



1=1 



Xe 



Vi 



N-l 



X g X, yi , 



yg + Xi 



VnVn-iXn-i 



xe-i ye-i + xe-i J 2<e<N-i ' ^Af-i(yAf-i + ^w-i) 
Change variables in the inner integral from x to y2.N so that the (7-integrand becomes simply g(x^ y). Recall 
that 7„ j = 6j + On- After matching up all the powers of the variables the integral above acquires the form 



n 



dx dy exp 



N-l 

E 



vi+i , yi+1 , Xi 



e=i 

, '^".j „~, -In, N-l 



yi 



XI 



N 



I ^ ye 



e=i 



N-l N 

X n ^f"'^^' ■ 11(2^7"'^ ^r"'"") -aii^yy 

j=i i=i 

That this agrees with the integral coming from the left-hand side of (5.3) is just a matter of substituting in 
the explicit formulas of the kernels. The proof of Lemma 5.1 is complete. □ 

5.2. Proof of Theorem 3.6. Part (i) follows from Proposition 3.4 and the theory of Markov functions (see, 
for example, [-VI]), as this was described in Proposition 3.3. Condition (i) therein is readily verified, while 
condition (ii) is the content of Corollary 3.5. 

Part (ii) follows from part (i) and the definition of the kernel {y^ , ■). 

Part (iii) follows from (3.26) setting f{z) = Y\^=iXi{z)~'^^ via the use of relations (3.14), (3.15), (3.22), 
(3.23): 

N , „ N 



1 



U'iLly^ "^'^e iy) ■^'^^tJi 



f nx,(z)-^'M,^(y,dz) ^ -i— / flx^izyO^-'^M^iy^dz) 



my) 
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In fact, if the A,; have non-zero real part then / is not bounded. However if we set fn{x) = /(x)l„-i<^<„ then 
the claimed formula follows from dominated convergence and the fact that /.j.^ K^{y, dz) \ f{z)\ is bounded. 

Part (iv), equation (3.27) follows by integrating equation (3.26) over fjL^ {y'^,dy) or using the intertwining 
relation (3.8) directly, as discussed in Remark 3.7. 

To deduce (3.28) we use the Plancherel formula (3.19). Let jj,^ iy'^,y) denote the density of {y^,dy). 
It exists because ^^{y^,dy) comes from composing kernels with densities. Multiply both sides of (3.27) 

/(o cx))" fiv)^^ {y)^^xiy)Y\iyi~^^Vi^ integrate over lMJ^ with respect to si^{\)d\ and use (3.19). The 
application of the Plancherel identity is valid since f{y)'ii^{y) is in L^((0, oo)^, Hi c^yj/j/i) for / continuous 
and compactly supported on (0, oo)^, and by the next lemma. 

Lemma 5.2. Let = 0i + 0j > be a solvable parameter matrix. Assume 6j — di < < 6i for all i,j > 1. 
Then for all y^ G (0, oo)^ and n > 1, the function 

IS m L2((0,oo)^,n,%/2;»). 

Proof. Iterating definition (3.11) and using (3.15) gives 



where (y°,y) is the density of the n-fold composition of kernels P^m{y^,dy^), . . . ,P^„](?/" ^,dy) from 
(3.6). Let 



ra{u,v) = T{a)~ (-) , -u, ue(0, oo), 

\V / V 

denote the transition density of a multiplicative random walk on (0, oo) with r^^(Q;)-distributed steps, and 
let 

» n 

i?„j(u°,u")= / TTr^. .(u*-\ u*)du"-i •••du^ 

denote the n-step transition density with parameters "fi,j from column j of our solvable matrix. By dropping 
the killing term from (3.6) and by an application of Jensen's inequality, 

N _ N 

{p:!{y",y)Y < lii^-Ay^^y^^y ^ n / i?n-i,,(y°,%)K„fe, 

.7 = 1 .7 = 1-^0 



J = l 3 = 

N. 



Put these together, noting that the integrals factor over (0, oo) 



i=i 

< '7 n / yT'^^Rn-Uj{y'j,yj){rj^,,{y], yj))^ dyjdjjj 



n / yf'Rn-h]{yj,y3)dyj 



^9 [y I "'(0,cx)) 



J=l 

=^w)^nnm.] < oo. 

The first equality above integrates away the variables yj. and the finitcness of the constant (7Ar.„(^["1 , 6*) 
depends on On > 0. The second equality uses the independence of random walk steps. The finitcness of the 
expectations is equivalent to 9j — 9i < 0. □ 
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5.3. Proof of Theorem 3.8. To prove Part (i) set y° = in (3.28). Asymptotic relation (20) in [42] 
gives that, as M — >■ oo, the ratio VPj^ (2/°'*^)/^'e' (2/°'*^) — ^ 1- Therefore, we only need to demonstrate that 
the limit M — > 00 can be passed inside the integral. This follows from dominated convergence since 



M/A'(j;0,A/) 



< 



*o^(2/°'*0 



for A S lR^ and the fact that the rest of the integrand is in L^{iR^ , s^r (A)dA). The latter follows from the 
fact that J(^Q f{y)'^e{y)'^-x{y)Y[i'^y'i/yi^ L'^{lR^ ,SN{X)dX), by the Planchcrel isomorphism and the 

fact that / is bounded and compactly supported, and nmnir(^m + K)/^{()i + ^m) G L^(tR^, sjv(A)ciA), 
for n > N. Indeed, using the asymptotics 



lim |r(.Ti +tX2)|e^l"^l|a;2| 



/27r, 



xi,x2 e 



it follows that 



N 



T{e,n + A,) 



sn{X) 



-A, I 



< 



-™E"iIa.|+^(w-i)E.=iIa.I 



which decays exponentially when n > N. 
Part (ii) follows from the fact that 

/ e-'y'^e{y)^-x{y)T\ 

"'(0,00)" 



dyi 

yt 



) e L^{iW\sN{X)dX). 



We can now repeat the argument used in the proof of Part (i), using the function f{y) = 6^*^^ (instead of 
a compactly supported /) and invoke the Whittaker integral identity (3.20). 

To prove part (iii) we show in Proposition 5.3 below a more general statement: the entire array z(n) 
converges in distribution, as M — >■ 00, to the one defined by the path configurations, noting that this 
statement can make sense only for the portion of the array constructed by time n. 

Let denote the probability distribution of the process z(-) when the initial state is ^(0) = z E Tn, and 
let denote expectation under . Let us also use the notation when the array starts empty, in which 
case at time n only the portion {zki{n) : l<k<N,l<£<kAn}oi the array has been defined. Recall 



y 



0,M 



'KKN 



with p from (3.30). 



Proposition 5.3. Let N,n> 1, and let f he a hounded continuous function of the (0, oo)-valued coordinates 
{zkiis) : 1 < s < n, I < k < N, I < i < k A s}. Then 



(5.6) 



lim 

M-)-oc 



(y"'^« , dz) E^[f{z{l), . . . , z(n))] = i?"[/(z(l), . . . , zin))]. 



Before turning to the proof of Proposition 5.3 we use it to derive part (iv) of Theorem 3.8. Let v'^\dz) 
denote the initial distribution {y^'^^ , dz) on arrays, and let / be a bounded continuous function on arrays. 
Then part (iv) follows from this calculation: 

E^[f{z{n))] = hm E'^'"'[f{z{n))] = lim i?"^" / {y{n),dz) f{z) 

= E^ f K^{y{n),dz)fiz)^ f {dy) f {y , dz) f {z) . 

The first and third equalities are instances of (5.6), the second is (3.25), and the last one is Proposition 
5.3 again because one consequence of limit (5.6), together with fi^ {y^'^'' ,dy) — )- fj.'^ [dy) from part (i), is 
that is the marginal distribution of y{n) under P** for n > N . The third equality above is justified 
by arguing that J (y, dz) f{z) is a continuous function of y, or equivalently, that y 1— >■ K^{y,dz) is 
a continuous mapping into the space of probability measures on arrays (in the usual weak topology of 
probability measures, generated by bounded continuous functions). This follows from the fact that, off the 
bottom row, (y, dz) has a density that is jointly continuous in (y, z). Pointwise convergence of densities 
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implies convergence of probability measures, a result known as ScliefFe's theorem. This completes the proof 
of part (iv). 

To prove part (v) of Theorem 3.8 we will use the alternative form of the Whittaker integral identity (3.21). 
By part (i) of the theorem, we have 



N N 



1 Tie, + e,„.) 

for any continuous, compactly supported function / on (0, oo)^. By (3.21) 

N , N N 

As above, the functions U /(y)*^(2/)*^!;,(y) Ui IT and 11^1=1 UZi n^+O^) arc both in L^lR^ , sjv(A)dA) 

so we have, by the Plancherel theorem, 

r N N „ N , 

lf{y)f^%{dy)=l[ llm+Lr' /^/(2/X(2/)e-^~ *f(y)n— ' 

•'^+ m=li=l •'^+ 1=1 

as required. 

Proof of Proposition 5.3. Figure 2 shows that {zke{s) : l<s<n, l<k<N,l<£<kAs} can be written 
as a function of {zm+iim) : < m < nAN} and rfl^'"'. Let (z(l), . . . , z(n)) = G{^{zm+i{rn))o<rn<nAN, d^^'"') 
represent this functional relationship defined by the row insertion procedure. The case of starting with 
an empty array is the one where each vector Zm+iiiri) = where e^*^"* = (1,0,..., 0) is the first k- 

dimensional unit vector. This can be seen from Figure 3. If we let P denote the probability distribution of 
the weight matrix d, the goal (5.6) can be re-expressed as 

1 



(Notation P(d((i[i'"l)) means that the matrix is the integration variable under the measure P.) On the 
left above the vectors {zm,+i{fn))(Xm<nAN sltc themselves functions of the initial values {zm{0))o<rn<nAN = 
{zm)o<m<nAN and thc weights as shown in Figure 2. Comparison of the ^' output in equations (2.1) 

and (2.2) shows that the mapping G is continuous as thc inputs z„i+i{m) — !■ e^^ Thus thc upshot is 

that we need to show the weak convergence {zm+i{m))o<m.<nAN ^ "^^)o<m<nAN as M oo. Since 

the limit is deterministic we can ignore the joint distribution and do this one coordinate at a time. So it 
suffices tofixO<m<A;<iV such that m < n A N and show that 

(5.8) 2;fc.m4-i(m) — > 5k,m+i in probability as Af — >■ oo 

when Zk^m+i(jn) has the probability distribution described by the left-hand side of (5.7) and 5k,m+i is the 
Kronecker delta. 
Write 

(5.9) Zk,m+l{m) = Vk^m[{zk,t[0))l<t<m+l,l<k<N,d^'^'"^^] 

to indicate the functional relationship from the inputs to the array element Zk^m+iiTn). Our goal (5.8) follows 
if we show that for any fixed dl^'™! £ (0, oo)"*^ and a bounded continuous test function /, 

(5.10) lim / iff(y°''*^dz)/(Ffe,,„[(zfc,,)i<,<,„+i.,<fc<jv,d[''"'l]) =/(4.n+i). 



Jim_ / 7^e^(y"^^^dz) / _ P(d(d[i'"l))/(G((z„+i(m))o<™<„A7V,d[''"l)) 

(5.7) 

P(d(d[i'"l))/(G((e(^-"))o<™<„A^,d[i'"l)). 



To understand the asymptotics of {y^'^^ ,dz) it is convenient to switch from multiplicative to additive 
variables. Define an array t ~ {tki}i<e<k<N by Zke = e**"^. Let Kly, dt) denote the distribution of the array 
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t when z has distribution K^{y,dz). For u E M.^ let 

W{u) = {t = {tki}i<i<k<N ■■ tN,t = Ui, 1 <i < N} 
be the set of arrays with bottom row u. Let 

N k-1 k N-1 k 

(5.11) Mt) - E ( E ^^-1/ - E ^^^^) -EE {e'^''-'^^''' + 6*'-+^.^+-*'=.^) . 
fe=i f=i f=i fc=i f=i 

Then, for a bounded continuous test function g, 

[ git) Kiy'-'', dt) = -1- / g{t) 6^"^ dt 

Above C{M) — /yj/j-gj e'^"'-*-'''"^"^^-^"^*'' dt is the norniahzation needed for a probabihty measure. We changed 
variables by shifting t to t — Mp where p = {pki)i<e<k<N is the array from (3.30) defined by pki = ^{k — 
1) - ^ + 1. Defining Se{t) = Fe{t) - ^"0(0 leads to 

Fg{t - Mp) = Sg{t- Mp) +Mt~ Mp) = Se{t) + e^'/^Mt)- 
Return to the right-hand side of (5.10) to rewrite as 

Xe^(y°'^^dz)/(Ffc.™[(zfc,,)i<f<„,+i,,<fc<jv,d['^"l]) 

= yi^(y°^*^d^)/(F,,™[(e*'=0l<^<m+l,^<fc<7v,d^''"l]) 
1 

C{M) 



(5.13) = -i- / /(Vi.,„,[(e*--^^^-)i<,<„,+M<..<^,d[i-™l]) e^-'W+^^'^^^W dt. 

lw{o) 

We claim that 

(5.14) line (5.13) converges to f{6k,m+i) as M — > 00. 

Limit (5.14) finishes the proof of Proposition 5.3. To establish it we prove the two lemmas below. 

Lemma 5.4. On the set W{0), the function J-q is strictly concave and has a unique maximum that 
satisfies ~ ^ f^'^ each 1 < k < N . 

Lemma 5.5. For each fixed 1 < m + 1 < k < N , m < n, dli-'"! G (0, 00)'"^, 

(5.15) lim Zk,m+i{m)= lim Ffc,„[(e*'=*-*^'"=Oi<^<™+i/<fc<iv, d''^"'] = <5fc,m+i. 

J\f— >oo M— >-oo ~ ~ ~ ~ 

Lemma 5.4 implies that the probability measure C(M)~^lvi/(o) (0 s'^"^*''"'''^"^^'^'''"*'' converges weakly to 
the pointmass at This together with (5.15) and the boundedness and continuity of / imply (5.14). We 
have proved Proposition 5.3 but it remains to prove the lemmas above. □ 

Proof of Lemma 5.4- This lemma comes from [4'2\ and [50]. We include the proof for the sake of completeness. 
The critical point equations -f^J^oit) = for I < i < k < N rearrange to 

e^'"' = , for 1< i < fc < . 

e-'''-'■^-'l{k>^>l} + e-*"*''^ 

The case k = 1 gives tu = (^21 +i22)/2. Multiplying equations together gives 

g2E gE^Ji'tfc-i..+E^^i'*k+i.i+i for 2 < fc < iV. 
From this follows fc~^ J2'i=i ^ki = J2i=i ^Ni ior 1 < k < N, and these all = by the W{0) condition. 
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Following [50, p. 136] we write in the following form. Consider the directed graph (V, 5) with vertex set 
V = {(fc, i) : 1 < i < fc < -/V} and where £ contains all possible edges ((fc+1, i), (fc, if) and ((fc, i), (fc+1, i + 1)). 
Edge a = (u(a),i;(a)) is directed form vertex u(a) to vertex v{a). Then 



gtii(a)-t„(a) 



j-o(t) = -E< 

and for vertices a;, y 
dtxdty \ 

Take a vector (Q!:c)2;gv such that at^f^.i) ~ (because the variables t^.i are not free to vary on W{Q)). Then 

is < unless a = 0. This gives strict concavity and the unique maximum. □ 

Proof of Lemma 5.5. First we take care of the case m = 0. This is read off directly from the initial values 
and = 0: Zfe,i(0) = e*k.i+M{i-(k+i)/2) _^ g^^^ as M ^ oo and t ^ 

For the rest of the proof m > 1. We turn to the matrix machinery developed in [41]. For that purpose 
we consider the row insertion procedure also in terms of ratio variables. Let rj = (rjg, . . . ,ri]y) denote the 
ratio variables associated with the vector <^ = (^£, . . . , ^^v): S,k = ViVt+i ' ' ' Vk for k = . . . ,N. Similarly 
^k ~ ViV'i+i ' ' 'V'k- Then the row insertion 

b b 

(5.16) ^ defined in Definition 2.1 is equivalently expressed as r] r]' . 

b' b' 

RecaU definition (2.13) of the TV x TV matrices Hm{v)- Then (5.16) is equivalent to [41, eqn. (2.23)-(2.25)] 

(5.17) H,{ri)H,{b) = Ht+^{b')H,{rj'). 

In the extreme case £ = N there is no b' left and the correct interpretation is HN+i{b') ~ I ~ the N x N 
identity matrix. 

As in the end of Section 2, define the ratio variables of the arrays by 

(5.18) rjuin-) = zu{n) and ■qhi{n) = Zki{n) / Zk-i/{n) ior 1 < i < k. 
Applying (5.17) to the upper left corner of Figure 2 gives 

i/i(77i(0))i?i(ai(l)) = i/2(a2(l))i?i(m(l))- 

Left multiply this identity by i^2('72(0)), ^^3(773(0)), . . . , right multiply by Hi{ai{2)), Hi{ai{3)), . . . , and 
apply (5.17) repeatedly on the right-hand side. This gives the following identity for all m > 1: 

m rn m m 

(5.19) ]^i?m+i_i(7?m+i_i(0)) • J]^ Hi{ai{j)) = J]^ Hm+2{am+2{i)) • ]^ffm+l-i('7m+l-i(n^)) 

1=0 j=l j = l j=0 

If 771 = — 1 then i?m+2(aTn+2(j)) = I and the first product on the right disappears. If 771 < A^ — 1 then 
apply [41, Thm. 2.4] to the lower right (A^ — 771 — 1) x (A — 771 — 1) block of the first product on the right. 
This gives vectors . . . such that = {pI, ■ ■ ■ ,Pn), 

m £-1 

n^'"+2(«-+20-)) =i[H,M,.,{p"^+'+^-^), 

J = l 3=0 

and £ = m A (A — rri — 1). Substituting this back into (5.19) gives 

m m £—1 m 

(5.20) llH,n+l-^{lU+l-^m ' J] = 11 H,n+e+l-j{p"'+'+'-') • [] (m)) • 
1=0 j=l j=0 i=0 
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Figure 5. The minor t™~^^ equals the sum of the weights of (m + l)-tuples of disjoint 
down- right paths (7^, . . . , 7™+^), where 7-' is a path from vertex j at the top to vertex 
k — m — 1 + j at the bottom. The thickset hne displays one admissible path 7^ from top 
vertex 2 to bottom vertex A: — m + 1. 



Let H denote the matrix on the left. On the right we have a descending sequence of subscripts (2m + 1) A 
A^, . . . , 1. We can appeal to [41, Prop. 1.6] to conclude that the vectors on the right-hand side are uniquely 
determined. In particular, 77m+i(m) = ("1))^=™+! is given by 

(5.21) Vni+i,,n+i{m) = , r]k,m+i{m) = J^~^ for m + 1 < fc < TV, 

■^ni+l '^k ^k-1 

where r] = detijl^'''. , -, , . i < 7, are minor determinants of the matrix H over rows 1,. . . ,i and columns 
j — i + I, . . . , j . Switching back to z via (5.18) gives 

(5.22) zk,m+i{m) = ^^ ioTm + l<k<N. 

This is the function Vk.m defined in (5.9). (Prop. 1.6 of [41] needs a hypothesis on the minors of H. This 
hypothesis can be checked from (5.23) below with the help of Figure 5.) 

We use a graphical representation to compute the minors rj, in the spirit of the Lindstrom-Gessel-Viennot 
method, following Sect. 1.1 of [41]. The matrix H is represented by an array of 2m + 1 right-adjusted rows, 
one row for each vector rj^+iiO), ■ ■ ■ , 771(0), d'^l, . . . , d'™' (see Figure 5). For 1 < i < m + 1, the vertices 
on row i are assigned weights ri^^2-i.m+2-i{0), ■ ■ ■ , ?77v,m+2-i(0), and for m + 2 < i < 2m + 1, the vertices 
on row i are assigned weights (ii-m-1,1, . . . iC^j-m-i.Af. Note that due to initial elements missing from the 
7y-vectors, the top vertex of column j is on row m — j + 2 for 1 < j < m. Combining (1.16) and (1.27) in 
[41] gives 

(5.23) r; = detiJf[;:^I^,^^.] = ^ wt{^\...,^^) 

(7^---,7*) 

where the sum ranges over i-tuples of disjoint paths 7^, ... ,7' such that 7'"" goes from vertex k at the top 
edge of the graph to vertex j ~ i + k a.t the bottom edge, and the weight wt{'-f^, . . . , 7*) of the i-tuple is the 
product of the weights on the vertices of the paths. 

Now we find the asymptotics of the minors in (5.22). We are proving (5.15) so the initial z^^-values are 
Zke = 6*'"*"^^'"='' = gtfcf+A/(«-(fe+i)/2)^ From this we get the initial ratio variables 

77« = zu = e*«+^^(^-i)/2, and rjue = Zki/zk-i,e = e*'^^-*'=-i,^-W2 fo^ k > £. 

Note in particular that on the top m 4- 1 rows of the array in Figure 5, all but the left edge weights decay 
as 

Consider first t™+^ for some fc > m 4- 1. One can check by induction on m that every (m + l)-tuple of 
paths (7^, . . . , 7™+!) from vertices (1, . . . , m 4- 1) on the top edge to vertices (fc — m, . . . , fc) on the bottom 
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Figure 6. An m-tuple of paths for the minor r™ 



edge contains at least m(m + l)/2 + 1 vertices with weight rj^i with £ < k. Since the t^i variables converge 
to a finite constant, up to a constant multiple 



«;t(7\...,7"+l) <Cl[ e*^(^-l)/2 . (g-Af/2)m(r„ + l)/2+l < Cg-A//2 



£=1 



There is a fixed finite number of these (m + f )-tuplcs in (5.23), and so r™^^ < Ce^^^/^ for fc > m + 1. 

Next we establish a lower bound for r™. The m rows of d- weights at the bottom of the array allow an 
m-tuple of paths that uses exactly m{m — l)/2 vertices with weight "qke with £ < k (Figure 6). This m-tuple 
gives a positive constant lower bound: t™ > c > 0. 

Combination of the first two bounds gives 



m+l 



(5.24) 



^:fc,m-i-l(m) 



< ^0 for m 1< fc < iV. 



It remains to consider the case fc = m + 1. Minor t^^I has a unique admissible (m + l)-tuple, namely 
m+l vertical paths. Consequently 



m+l 
'm+l 



n 



Vke ■ 



n "^^j 



n ^^^-^ 



l<£<fe<m+l 



l<i<rn 
l<j<m+l 



l<i<m 
l<i<m+l 



For T^^i the m-tuple in Figure 6 has weight 



l<e<k<m l<i<m 
l<j<m+l 



n 



l<2<m 
l<:/<m-|-l 



This m-tuple is minimal in its use of weights rjki with £ < k. Other admissible m-tuples for r™^^ necessarily 
use more of such weights and consequently pick up more e~^^/^ factors. From all this 



Zm+l,m+l{lTl) 



(5.25) 



'm+l 



'm+l 



n l<j<m di.j 
l<j<m+l 



^E.-Li . n i<.,<m + 0(e-*V2) 

l<i<m+l 



^m + l ^0 



1. 



(5.24) and (5.25) together verify (5.15) and complete the proof of Lemma 5.5. 
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